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FIG. 4 Plan view of input buoyancy structure and faufts (a), model output of 
surface velocity (red arrows) and <6), topography and long-wavelengtn 
{> 1.000 km) geoid (c) for the model with variable strikes of slabs, and 
surface velocity and <d) for the model excluding the faufts. Contour 
interval for geotd is 10 m and dashed lines are negative values. The colour 
W scales are the same as m Fig. 2. The rengtn. width and thickness of the 
model boxes (in km) are 6,000, 7,500 and 1.500, respectively; which are 
covered by 54 x 54 x 16 elements. 



faults and slabs with variable strikes (Fig. 4a) shows that surface 
velocities remain plate-like (Fig. 46) with P = 0.87 and 
R T ? = 0.40 {R T;? = 0.42 for a perfect plate). Surface velocities 
in R are predominantly parallel to the vertical faults and show 
oblique subduction where the dipping faults are not parallel to the 
spreading centre (Fig. 46). Trench topography is evident ocean- 
ward from the dipping faults and there is 40 m long-wavelength 
geoid high over slabs (Fig. 4c). Qualitatively, both trench/fore- 
bulge topography and the geoid are essentially identical to 
observed features in the western Pacific (for example, approxi- 
mately 4 km deep and 150km wide trenches, and 30m long- 
wavelength geoid high). The reduced P results form the velocity 
variations in the spreading perpendicular direction near the short 
vertical fault (Fig. 46) where slabs are younger; a more realistic 
temperature-dependent rheology may lead to a higher plateness. 
When the faults are excluded from the model, surface velocities 
vary gradually in both the spreading parallel and perpendicular 
directions with P decreased to 0.75 (Fig. Ad). However, R T „ only 
decreases slightly to 0.39, suggesting that a large fraction of 
toroidal power is not a sufficient criterion for characterizing the 
existence of plates. 

These models show that the characteristics of plate motion can 
be explained with realistic mantle buoyancy, weak faults and a 
power-law rheology (n = 3). Besides plate kinematics, the models 
also simultaneously predict realistic trench/fore-bulge topography 
and the long-wavelength geoid (Fig. 46 and c). In the past, 
different physical processes were used to explain these observed 
features (for example, elastic plate bending models for trench/ 
fore-bulge topography and laterally homogenous viscous flow 
models for the geoid). By simultaneously matching observed 
trench topography, plate kinematics and the long-wavelength 
geoid, our models will provide more realistic constraints on 
materials properties and slab dynamics. Weak transform faults 
enhance the motion along the faults. Weak transform faults and 
their strong adjacent regions cause transform faults to guide plate 
motion (Fig. 2e and h), consistent with inferences based on 
plate kinematics 13 . The guiding influence of transform faults on 
plate motion and the decoupling of thrust faults result in oblique 
subduction where the strike of thrust faults is non -perpendicular 

NATURE • VOL 383 • 19 SEPTEMBER 1996 



to rhe transform faults (Figs 2h and 46), consistent with observa- 
tions 1 - 1 *. Non-zero normal velocities on faults indicate migrating 
plate margins. Ultimately, with the formulation presented here 
and with multiple plates coupled with the solution of energy 
equation, we should be able to predict the evolution of plate 
geometry as accomplished in rwo dimensions' 1 . □ 
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Asthma now affects one child in seven in the United Kingdom 1 . 
Most cases (95%) of childhood asthma are associated with atopy, 
the immunoglobulin E (IgE)-mediated familial syndrome of 
allergic asthma, eczema and rhinitis. Segregation analysts has 
consistently suggested the presence of major genes influencing 
atopy and IgE levels^, with the expectation that these genes may 
be identified by positional cloning or the examination of candi- 
date genes. Here we report the results of a genome-wide search for 
linkage to one qualitative and four quantitative traits associated 
with allergic (atopic) asthma. We have identified six potential 
linkages (P < 0 001), five of which are to quantitative traits. 
Monte Carlo simulations show that 1.6 false-positive linkages 
at this level of significance would be expected from the data. One 
linkage, to chromosome UqI3, has been established previously 5 . 
Three of the new loci show evidence of linkage to a second panel of 
families, in which maternal effects and pleiotropy of linked 
phenotypes are seen. The results demonstrate the extent and 
the complexity of the genetic predisposition to asthma. 

Asthma is usually recognized epidemiologically by standard 
symptom questionnaires or by physician diagnosis 6 . The non- 
specific bronchial hyper-responsiveness that accompanies 
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asthma is measured by the dose 
response of airflow to a broncho- 
constrictor such as histamine or 
methacholine. Atopy is detected 
by skin prick tests, by measure- 
ment of specific serum IgE titres 
against allergens using radio- 
allergosorhent test (RAST) or 
enzyme-linked immunosorbent 
assay (ELISA) techniques, or by 
quantifying the total serum IgE. 
Atopic disease is accompanied by 
blood eosinophilic and eosinophils 
are prominent in asthmatic 
airways. 

The prevalence of atopy in 
Western populations is between 
40 and 50% 7 ~". The recruitment 
of families more than one 
asthmatic member results in 
samples in which 10% or more 
subjects are atopic, making it 
much mAffe difficult to detect 
linkage 10 For this reason, the 80 
families in the genome screen 
(from Busselton in Western 
Australia) were subseiected 
from a population sample to 
include non-atopic members. 
We also reasoned that detection 
of linkage to quantitative traits 
would be enhanced by including 
subjects in normal and abnormal 
ranges of the trait distributions. 

We looked for linkage to four 
quantitative parameters: the total 
scrum IgE. the skin test index 

(STI), the peripheral blood eosinophil count, and bronchial 
responsiveness to methacholine (slope)". A RAST index was 
also calculated, but gave no additional information to the STI 
and so was not included in the analysis described below. 

Differing indices of atopy may be increased in the same 
family 12 -", and RAST and skin-test responses reach a peak later 
in childhood than the total serum IgE, and decline at a slower rate 
thereafter". To .account for this heterogeneity of phenotype 
(pleiotropy), the categorical trait of 'atopy' was used, based on a 
combination of STI, RAST index, and total serum IgE. Linkage to 
quantitative traits was tested by the Haseman-Elston sib-pair 
technique 14 , and to atopy by affected sib-pair methods. 

The 80 families surveyed contained a total of 203 offspring 
forming 172 sib-pairs. The mean age of the children was 12.6 ± 1.3 
(s.e.) years, their geometric mean IgE was 55.7 db 1.1 kUI~\ and 
their mean STI was 4.0 ±0.41 mm; 12% of the children were 
asthmatic, and 25% admitted to wheeze; 52% were atopic, 36% 
were non-atopic, and 12% had an intermediate phenotype as 
defined; the latter were classified as unknown. The sib-pairs 
subdivided into 44 pairs in which both siblings were atopic, 87 
pairs in which one sibling was affected, and 41 in which both were 
unaffected or of intermediate phenotype. 

The families were initially screened with 269 markers (253 
autosomal and 16 X-linked). A genome map was constructed 
from these data, and the order and distances between the markers 
were compatible with previously published values 1516 . Based on 
these families, the markers had an average heterozygosity of 75%, 
and spanned approximately 3, lOOcentimorgans (cM). After initial 
statistical analysis, none of the X chromosome markers showed 
evidence of linkage to the five phenotype traits, and the remainder 
of this analysis therefore deals exclusively with autosomal markers. 



TABLE 1 Positive results of genome search (empmcal P values) 
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Only regions showing linkage to one or more phenotypes with P < 0 001 are shown. The recombination 
fraction ((?) is to the subsequent marker. 



somes 4, 6, 7. 11. 13 and 16, respectively (Table 1). For the regions 
on chromosomes 4, 7 and 16, the test statistic met the more 
stringent criterion of P < 0.0005. and for two regions, on chromo- 
somes 6 and \ \,P was <0,0001. Removal of sib-pair differences 
that were greater than 3 standard deviations from the mean did 
not significantly alter the P values for these loci. 

The preliminary* identification of regions of potential linkage 
used correlated phenotypes with markers spanning most of the 
genome. Because of multiple testing, the true significance values 
will be greater than the normative values given in Table 1. To 
judge the global significance of the results, a Monte Carlo 
procedure was applied to calculate the number of regions that 
would be expected to reach these levels of significance in the 
absence of any linkages (that is, the number of false -positive 
linkages). The procedure has advantages over traditional 
approaches to correcting nominal significance levels because it 
allows for the in format iveness and distance between markers that 
were actually characterized in the genome-wide search. The 
expected number of regions of false-positive linkage was 1.6 for 
P < 0.001 (compared with 6 observed linkages), 0.8 for 
P< 0.0005 (5 observed linkages) and 0.2 for P < 0.0001 (2 
observed linkages), so the observed number of regions meeting 
these statistical criteria for linkage was much larger than would be 
expected by chance. In 35,000 simulations, P < 0.001 was 
observed in less than 2% of replicates, indicating that the prob- 
ability of seeing six or more such regions without a true linkage was 
<0.005. 

Additional markers were studied on the six chromosomes that 
showed evidence for linkage (Table I). Although the evidence in 
favour of linkage was not increased, in all instances at least two of 
the markers from each region had significance of P < 0.05. 
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the p-chain of the high-affinity receptor for IgE 
(FceRI-P), also showed linkage to the STI, but this 
was not observed on chromosome 16. The regions on 
chromosome 4 and 7 were linked to bronchial respon- 
siveness, whereas the region on chromosome 6 (near 
the class I genes of the major histocompatibility com- 
plex (MHQ) was linked to eosinophil counts. Weaker 
evidence of linkage {P < 0.01) was found with other 
phcnotypes for the markers on chromosomes 6 and 7 
(Tabic 1). 

Markers from chromosome 13, around D13S153, 
showed evidence of linkage to the atopy phenotype in 
affected sib-pair analyses (Table 1). The sharing of 
maternal and paternal alleles identified by descent 
(IBD) was 46,1 and 20,0. 

The markers showing P < 0.001 for linkage were 
tested for replication in an additional group (panel B) 
of 77 nuclear and extended families, recruited from clinics in the 
United Kingdom, who had previously been used to map atopy on 
chromosome llql3\ These families had been characterized by 
the same protocols as the Busselton families, except that bronchial 
responsiveness and eosinophil counts had not been measured. 
They contained 215 offspring (268 sib-pairs), of which 61% were 
atopic ancT 56% asthmatic, the high values reflecting their origin 
through clinics. There were 121 atopic-affected sib-pairs, and 70 
asthmatic sib-pairs. 

Linkage of asthma to FCERB (P = 0.003) and to D16S289 
(P = 0.03) was seen in the panel B families (Table 2). Linkage 
with atopy was found to D13S153 (P = 0.003). The sib-pair 
sharing of D13S153 alleles in these subjects was 65 IBD = 1, 
36 IBD = 0 (P - 0.003). The sharing for the combined data was 
1 1 1 .IBD = 1, 56 IBD = 0. The combined Hodge-weighted statis- 
tic was 4.158 (P = 0.000016). The empirical P value for linkage to 
D13S153, calculated after 500,000 estimations, was <0.00001. 

Many studies have shown that maternal atopy and asthma are 
more likely to be transmitted to children than is paternal 
atopy' 7 :o We have not tested for maternal linkage in the 
genome screen families because of the difficulties induced by a 
further set of multiple comparisons. However, maternal linkage of 
atopy to chromosome 1 Iq markers has previously been demon- 
strated in the panel B families 510 , and so potentially linked 
markers identified in the genome screen were tested for parent 
of origin effects (Table 2). 

D4S426. FCERB and D16S289 showed significant differences 
in linkage between maternal and paternal alleles (Table 2). 
Linkage to maternal meioses was seen between D4S426 and 
atopy (P < 0.05) and the total serum IgE (P < 0.001), and 
between D16S289 and atopy (P < 0.01) and asthma 
(P < 0.001). As already observed, FCERB showed strong mater- 
nal linkage to atopy 5 10 , but a strong linkage to asthma was also 
seen (P < 0.00001) which has not been previously noted. The 
finding of maternal effects at several loci favours immunological 
interactions between mother and child, rather than genetic 
imprinting or anticipation, as a cause of these phenomena. 

ft is not clear if linkage of different traits to different regions of 
the genome indicates that particular genes affect specific traits. 
However, a previously established region of genetic linkage, to 
chromosome llql3, may be used as an examplar 1321 * 24 . In the 
genome screen, the 1 lq!3 marker FCERB was linked to the IgE 
and the STI. In other studies, FCERB or neighbouring markers 
have been linked to the composite phenotype of atopy 13 - 1 - and to 
bronchial responsiveness apparently independently of atopy 24 . 
Similarly, D4S426 and D16S289 showed linkage to different 
phenotypes in genome screen and the panel B families. Although 
D13S153 showed linkage to the same phenotype of atopy in both 
Australian and British families, linkage between the total serum 
IgE and a polymorphism in the esterase D protein, near D13S153, 
hasjDreviously been suggested, with a lod score of 2.07 (ref. 25). 



TABLE 2 Testing for linkage in panel B families: combined and maternal meioses 
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tropy may be attributable to environmental events, such as 
exposure to allergens, and to the age and selection of subjects. 
The relative distribution of quantitative traits in normal and 
abnormal ranges may also have differentially influenced the ability 
to detect linkage. 

The cytokine gene cluster on chromosome 5q31 has previously 
been reported to be linked to serum IgE in Amish families 2 *, and 
the result confirmed in Dutch asthmatic families : \ In our study we 
did not see linkage to either the serum IgE or to bronchial 
responsiveness with D5S393. D5S210 or D5S410 from the cyto- 
kine gene cluster, although we are currently typing multiple 
markers within this 20-cM region. 

In summary, both the Monte Carlo simulations and the replica- 
tion of some positive results in a second set of subjects suggest that 
we have identified regions of true linkage. Replication of linkage 
to complex traits is to be expected to be difficult 2 ", and testing for 
replication in other data sets and with a wider range of markers is 
still required for the linkages to chromosomes 6 and 7. It is 
unlikely that all important loci predisposing to atopy and 
asthma have been mapped, and larger studies and the pooling 
of data from several centres are required. □ 



Methods 

Subjects. The pnmary genome screen involved 364 subjects m 80 nuclear 
families subselected from a population sample of 230 families from Busselton 
in Western Australia* 29 . Clinical testing took place over three months in winter, 
to control for pollen exposure. Young families were recruited to reduce the effects 
of age and smoking on the IgE levels and lung function. Families in the genome 
screen included both atopic and non-atopic members, and sibships of three or 
greater were not exclusively atopic or non-atopic. 

Phenotypes. Skin tests to house-dust mite (HDM) and mixed grass pollen 
{minus the response of negative controls), specific IgE titres to HDM and 
timothy grass, and the total serum IgE were measured 29 . A "skin test index* (STI) 
was calculated as the sum of the skin-prick test results to HDM and grass mix; 
95% of individuals in this population who were atopic reacted to either H0M or 
grass pollen, or both). Bronchial responsfveness to (up to 12 nmol) methacho* 
line was measured 29 . The slope of the dose- response curve was calculated: 
(pre-dose forced expiratory volume in one second (FEV1) - lastFEVl)/cumula- 
tive dose of methacholine. A constant of 0.01 was added to each measurement, 
to allow loge transformation when the slope was <0. Eosinophils in peripheral 
blood were Coulter-counted and the values log, -transformed before analysis. 
Atopy was defined as STI > 5 mm, or a RAST score to HDM and timothy grass 
>2, or a total serum IgE > the 7th decile of the age-corrected population, 
'Normal' was defined as an STI of 0 and a RAST index of 0. and a total serum 
IgE < the 7th decile of the age- and sex-matched population. Intermediate 
phenotypes were classified as unknown. The subjects were administered a 
modified Bntish MRC questionnaire' 9 . Asthma was defined as a positive answer 
to the questions "Have you ever had an attack of asthma?" and "If yes. has this 
happened on more than one occasion?" 

Genotypes. Genomic ONA was isolated from penpheral blood leucocytes by 
phenol/chloroform extraction. We typed 22 dinculeotide repeat markers by 
radioactive techniques 15 (marker names italicized in Table 1). For 274 markers. 
50 ng of each DNA sample was geno typed by fluorescence-based semi-auto- 
matic methods 15 . Allele sizes were determined using the GENOTYPER soft- 



LETTERS TO NATURE 



man exact alleles were determined. Non-mendelian marker data were flagged by 
the programs GAS and UNKNOWN. Potentially incorrect genotypes were re- 
examined and retested if necessary. Two-potnt lod scores between markers and 
marker order were compared with published maps 1516 . 

Statistical analysis. Sib- pairs were included in tne linkage analysis only when 
marker genotypes were known for both parents. Linkage was sought to atopy by 
affected sib -pair analysis. Sibships containing more than one pair were 
weighted 30 . Linkage to quantitative trarts was calculated by the method of 
Haseman and Elston*\ in wnich the probabilities that sib- pairs have two, one 
or no alleles identical by descent are regressed o the square of the sib-pair trait 
differences 14 . Linkage was examined further in the region of markers where P 
was <0.001. For Table 1. linkage was presented with P < 0.01 or P < 0.05 if a 
marker is < 15 cM from a marker with a P < 0.001. The P values for linkage with 
the quantitative trait loct were assessed from the regression analysis of vanance 
(ANOVA), which assumes independence of the sib-pairs. In regions of potential 
linkage (P < 0.001). the P values were validated by empirical calculations in 
which this assumption is not required. The empirical statistics shown in Table 1 
were obtained by simulating the segregation of markers to offspring conditional 
on the observed parental genotypes without linkage to any of the traits. The 
emptncal P values are the proportion of instances in which the linkage statistic 
equalled or exceeded its observed value in 500.000 simulated replicates. The 
frequency of false- positive linkages was assessed by a Monte Carlo simulation, 
in wnich replicate sets of observations were simulated under the absence of 
linkage w)th mendelian segregation of the offspnng genotypes conditional on the 
observed parental genotypes. Linkage between markers was incorporated with 
recombination distances equal to those estimated from the data map obtained 
from this study^with offspring simulated under the assumption of lack of 
interference. The linkage statistics were calculated for each marker locus with 
the offspnng phenotype values as observed. The number of regions in which the 
linkage statistic for one or more phenotypes exceeded the assigned critical value 
was calculated over 35.000 replicates. Markers were assumed to detect 
different regions of linkage if all of the P values were non- significant over a 
region of >15cM between two markers with significant linkage. 

Replication. Panel B families were recruited from clinics in the United 
Kingdom 5 . Each family contained at least one asthmatic identified as a proband. 
Phenotyping was as for the Bussefton families, except that eosinophil counts had 
not been performed and bronchial responsiveness had not been measured. The 
families were tested with markers that showed potential linkage (P < 0.001) to 
any quantitative or qualitative trart in the Bussefton families. Genotyping and 
linkage analysis were as for the Bussefton families. In addition, loci were tested 
for the presence of significant differences in linkage between maternally and 
paternally derived alleles; if a significant difference was found, the significance of 
linkage to maternal alleles was estimated. 
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Anterior -posterior (A-P) patterning is of fundamental impor- 
tance throughout vertebrate embryonic development. Murine 
members of the trithorox (trx) and Poly comb group (Pc-G) of 
genes regulate A-P patterning of segmented axial structures 1 " 4 , 
demonstrating conserved upstream regulation of homeotic path- 
ways between Drosophiln and mouse. Here we report the posi- 
tional cloning of a classical mouse mutation, eed (for embryonic 
ectoderm development), which is the highly conserved homologue of 
the Drosophiln Pc-G gene esc (for extra sex combs), a long- term 
repressor of homeotic genes 5 . Mutants homozygous for a null 
allele of eed display disrupted A-P patterning of the primitive 
streak during gastrulation. Mutant embryos lack a node, noto- 
chord and somites, and there is no neural induction 4 . In contrast 
to absence of anterior structures, extra-embryonic mesoderm is 
abundant. Mice carrying a hypomorphic eed mutation exhibit 
posterior transformations along the axial skeleton. These results 
indicate thai eed is required globally for A-P patterning through- 
out embryogenesis. 

The recessive, embryonic-lethal mutation eed maps to the 
albino-deletion complex on chromosome 7 (ref. 7). Previous 
genetic and molecular analyses positioned eed within a 150- 
kiiobase (kb) interval between the c 3H and the c* L deletion 
breakpoints 7 (Fig, la). Here we have isolated eed from a 
CpG island (that is, a cluster of methylation-sensitive rare- 
cutter sites) located approximately 40 kb proximal to the c 4L 
deletion breakpoint 7 . Molecular access was obtained by 
chromosomal walking from the proximal end of yeast artificial 
chromosome (YAC) D6E7 (Fig. la). A 1.1-kbSafl fragment from 
the CpG island was used to isolate overlapping clones from a 
mouse embryonic day 7.5 (E7.5) cDNA library generating a 1.95-kb 
cDNA (Fig. la). The candidate cDNA was validated as the eed 
gene by analysing two non-complementing ethyl-W-nitrosurea 
(ENU)-induced alleles (l7Rn5 tiS4SB and l7Rn5' 989SB ) at the eed 
locus 8 * 9 . The l7Rn5 3}54SB and l7Rn5 ,Q89SB mutations represent 
null 10 and hypomorphic alleles 9 , respectively. DNA sequence 
analysis of reverse transcription-polymerase chain reaction 
(RT-PCR) products from homozygous E8.5 l7Rn5 m4SB embryos 
showed a transition of thymine at position 1040 to cytosine 
(T 1040 — C) (Fig. la). This base exchange was absent in an 
ENU-induced allele, l7Rn5 4234SB , which fww-complements the 
l7Rn5 3354SB allele 9 . Polymorphism can be excluded because both 
mutations were generated on the same BALB/cRl background 8 . 
As a second co-isogenic control, RT-PCR products from E12.5 
homozygous l7Rn5 r98gSB embryos also showed the expected wild- 
type T 1040 . The sequence alteration in l7Rn5 imsB was identified as 
a T 1031 — A transversion which was absent in both l7Rn5 4234SB and 
l7Rn5 m4SB (Fig. la). 

Conceptual translation revealed an open reading frame of 441 
amino acids, containing five tandemly repeated WD motifs (Fig. 
la). A potential PEST sequence" is defined by residues 22 to 52 of 
Eed protein (Fig. la). Two discrete clusters of basic amino acids 
encompassing a nine-residue spacer region (RK-X^-KXXKX, 
residues 65 to 80; Fig, la) are reminiscent of a bipartite nuclear 
targeting signal 12 . Eed does not display any canonical DNA- 
binding mom". Three putative Asn-linked glycosylation sites 
(NXS/T) 13 at amino-acid positions 157, 283 and 349 (Fig. la) 
are comoatible with a potential non-nuclear isoform of Eed 



